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Abstract 

Background: In plants, a large family of calmodulin (CaM) and CaM-like (CML) proteins transduce the increase in 
cytosolic Ca 2+ concentrations by binding to and altering the activities of target proteins, and thereby affecting the 
physiological responses to a vast array of stimuli. Here, transcript expression analysis of Cam and CML gene family 
members in rice {Oryza sativa L.) was extensively examined. 

Results: Cam and CML genes in rice exhibited differential expression patterns in tissues/organs. Under osmotic 
stress and salt stress, expression of OsCaml-1, OsCML4, 5, 8, and 7 7 was induced with different kinetics and 
magnitude. OsCML4 and 8 mRNA levels significantly increased by 3 h after treatment and remained elevated for at 
least 24 h while expression of OsCaml-1, OsCML5 and 7 7 was up-regulated as early as 1-3 h before rapidly 
returning to normal levels. Several c/'s-acting elements in response to abiotic stresses, including DREs (important 
promoter elements responsive to drought, high salt, and cold stress), were detected in the 5' upstream regions of 
these genes. The observed induction of the GUS activity of transgenic rice plants via the OsCaml-1 promoter 
appeared to be biphasic and dependent on the severity of salt stress. 

Conclusions: Large OsCam and OsCML gene family members likely play differential roles as signal transducers in 
regulating various developmental processes and represent important nodes in the signal transduction and 
transcriptional regulation networks in abiotic stresss responses mediated by the complex Ca 2+ signals in plants, 
which are rich in both spatial and temporal information. 
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Background task performed by various Ca + -modulated proteins. For 
Transient changes in the cytosolic Ca 2+ concentration the majority of these proteins, the Ca 2+ -binding sites are 
([Ca 2+ ] cyt ) of a different magnitude and specialized charac- composed of a characteristic helix-loop-helix motif called 
ter due to the activities of Ca 2+ -ATPases and Ca 2+ -channels an EF-hand [3], which binds Ca 2+ with high affinity, re- 
in the cellular membrane of eukaryotes are utilized as a suiting in conformation changes that modulate their acti- 
second messenger in generating physiological responses vity or their ability to interact with other proteins, 
to extracellular stimuli. The use of Ca 2+ signals has been Plants possess a complex network of Ca 2+ signal trans- 
implicated in generating responses to a wide variety of duction in mediating responses to various biotic and 
environmental changes in plants [1]. The changes in abiotic environmental stimuli. In plants, genes encoding 
[Ca 2+ ] cyt are not only transient, but also vary spatially EF-hand containing proteins have been extensively anno- 
and temporally with different organelles acting as dis- tated in Arabidopsis thaliana (L.) Heynh [4] and rice 
tinct compartments [2], therefore, the diverse array of {Oryza sativa L.) [5]. Three groups of EF-hand-containing 
different changes in the [Ca 2+ ] cyt must be discriminated Ca 2+ sensor proteins have been identified, which include 
so as to elicit the correct subsequent cellular response, a Ca 2+ -dependent protein kinase (CPK), calcineurin B-like 

protein (CBL), and calmodulin (CaM) [6]. These proteins 
can be considered as sensor responders or sensor relays 
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as they possess kinase activity either in their mole- 
cules or in their high affinity interacting partners as a 
responder function. By contrast, CaMs and CMLs, which 
have no other identifiable functional domain other than 
EF-hand motifs, are considered sensor relays. A large fa- 
mily of Cam and Cam-like (CML) genes has been exten- 
sively identified from the two model plants, A thaliana 
[8] and O. sativa [5]. Although the existence of these pro- 
teins in a single plant species is believed to be important 
to correctly perceive and discriminate the Ca 2+ signals 
from different stimuli, and thus aid in eliciting the correct 
subsequent response, the molecular mechanisms and 
physiological significance of most of these proteins have 
not been established. Nevertheless, accumulating evidence 
suggests that each of the different Cam and CML genes 
may have distinct and significant functions [9]. 

Currently, the information on Cam and CML genes in 
O. sativa is quite limited despite the fact it is a model 
plant for monocot and especially cereal plants. Regar- 
ding expression studies, reports on only a few members 
of the OsCam and OsCML gene family, are available 
[10,11]. Here, the transcript expression levels of several 
members of the OsCam and OsCML gene family were 
examined in O. sativa that had experienced abiotic stres- 
ses. Their 5' upstream sequences were examined to 
identify putative stress-responsive ds-acting elements. 
The putative promoter of the OsCam 1-1 gene was fur- 
ther examined using three transgenic rice lines. Based 
on their expression patterns, the OsCam and OsCML 
genes that are likely to represent important nodes in the 
signal transduction and transcriptional regulation net- 
works of Ca 2+ signaling in response to different abiotic 
stresses will be revealed, and subject to confirmation, 
this is likely facilitate the efforts in characterizing the 
mechanisms of abiotic stress Ca 2+ signalling in plants. 

Results 

Differential expression of OsCam and OsCML gene family 
members in various organs 

The rice (Oryza sativa L.) genome has been reported to 
encode for a large family of five Cam genes and 32 genes 
encoding CML calcium sensor proteins [5]. To deter- 
mine which of these OsCam and OsCML genes may 
have tissue- and/or stress- specific expression patterns, 
which could suggest their functions, their expression 
levels in different rice tissues and organs were investi- 
gated using publicly available microarray data of global 
gene expression from the experiment GSE 6893 from 
the O.sativa IR64 cultivar [12]. The available expression 
data of all OsCam and OsCML genes plus OsEFla were 
retrieved and the results are shown as heat map in 
Figure 1A, and are grouped according to the previously 
reported clades found by phylogenetic analysis of their 
amino acid sequences. It is clearly noticeable that all the 



OsCam genes (group 1) are highly expressed in almost 
all of the organs examined. In contrast, several highly 
conserved OsCML (groups 2-5) genes tend to be highly 
expressed in some organs/tissues. Interestingly, one 
member of each subgroup in group 6 is found expressed 
at a higher level compared with the other members in 
their subgroup (OsCML30 in group 6a; OsCML27 in 
group 6b; OsCML29 in group 6c; OsCML31 in group 6e) 
except group 6C, in which all members are highly 
expressed in most organs examined. 

Transcript expression levels of ten genes, one Cam 
(OsCaml-1) and nine highly conserved CML genes 
(OsCMIi, 3, 4 and 5 in group 2, OsCML8 and 13 in 
group 3, OsCML7 and 9 in group 4, and OsCMLll in 
group 5) were closely examined in different rice organs 
and tissues (Figure IB). Our preliminary results indi- 
cated that OsCaml-1 gene is the only OsCam gene 
which expression is induced by osmotic stress and salt 
stress, therefore, only OsCaml-1 among the five OsCam 
genes was selected for this study. Based on the micro- 
array data analysis, ubiquitous transcript expression at 
relatively high levels of OsCaml-1, OsCML7 and 13 was 
observed. The other three genes, OsCML i, 4, and 8 had 
similar expression patterns with higher levels in the leaf 
and the root and lower levels in the young inflorescence 
(Figure 1A). Their expression levels were then increasing 
upon maturation of the inflorescence until reaching 
similar levels to those of the leaf and the root in the 
early developing seed. Expression levels of OsCML4, and 
8 were maintained in the seed until dormancy and desic- 
cation while that of OsCML 1 later decreased upon seed 
maturation (Figure IB). For the other four genes: 
OsCML3, 5, 9, and 11, overall, they exhibited relatively 
lower expression levels compared to the aforementioned 
genes. 

To examine whether these genes are expressed in the 
three-week old Khao Dawk Mali 105 (KDML105) rice cul- 
tivar seedlings and to optimize the PCR conditions, RT- 
PCR was performed using the oligonucleotide primers 
specific for each gene (Table 1). Following separation of 
the PCR products in agarose-TBE gels and visualization 
by uv-transillumination after ethidium-bromide staining, 
bands of the expected size based on the OsCaml-1 and 
OsCML sequences (Table 1) were specifically detected 
(data not shown). Thus, transcripts of OsCaml-1 and all 
nine OsCML genes examined are expressed in leaves of 
the three-week old KDML105' rice seedlings, and can be 
specifically amplified by RT-PCR under these conditions. 

The transcript expression levels of these genes in the 
leaf blades/sheaths of the three-week old KDML105' rice 
seedlings, was evaluated by rt-RT-PCR using the trans- 
cript expression level of OsCMLS as reference (Figure 2), 
OsCaml-1 was found to have the highest expression 
level (~ 150-fold higher than OsCMLS) in the leaf tissue 
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Figure 1 OsCam and OsCML gene expression in various organs and developmental stages based on DNA microarray. (A) Expression 
profiles of the Com and CML genes in the 0. sotivo IR64 cultivar shown as a heat map. The groupings designated as 1-6 based on the 
phylogenetic analysis of CaM and CML proteins are indicated on the left and the colour scale is shown on the right. (B) and (C) Expression 
profiles of OsCaml-1 and the nine CML genes shown as bar graphs. Data in (B) and (C) are shown as the mean ± 1 SD, and are derived from 
three independent replicates. 



of the three-week old rice seedlings, followed by 
OsCML4, 8 (~120-fold) (Figure 2 A). OsCMU, 7, and 13 
also exhibited a relatively high expression level (>30-fold), 
while OsCML3, 5, 9, and 11 expressed at a relatively lower 
level 

The transcript expression levels of OsCam 1-1 and all 
nine OsCML genes assayed in the leaf blades/sheaths, 
roots, flowers, and seeds of the KDML105 rice cultivar 
during the grain filling period were determined and 
compared with their respective levels in the leaf tissue 
of the three-week old rice seedlings (Figure 2B). The 
expression profiles of these genes may be divided in- 
to three groups when only large differences (> 5-fold) 
among organs were considered. The first group, which 
has the highest expression levels in the root includes 
OsCML4, and 5. The second group with the highest ex- 
pression levels in the flower and the seed were OsCML3, 
and 11. Finally, OsCaml-1 and the five other OsCML 
genes (OsCML 7, 7, 8, 9, and 13) had more or less similar 
expression levels among these organs and comprised the 
third group. 



Expression of OsCam and OsCML genes under abiotic 
stress 

To investigate whether OsCaml-1 and the nine selected 
OsCML genes may possibly be involved in mediating 
responses to abiotic stress, the publicly available rice 
microarray data was retrieved and examined to acquire 
data for the IR64 rice cultivar, whilst the transcript and 
expression levels in the 'KDML105' rice were evaluated 
after different abiotic stresses by rt-RT-PCR. 

The results of the DNA microarray data of OsCaml-1 
and the nine OsCML genes in the 7-day-old 1R64' rice 
seedlings grown under 3 h of dehydration (drought) or 
salt stress (200 mM NaCl) (GSE 6901) [12] retrieved via 
the rice oligonucleotide array database are shown as a 
heat map in Figure 3. The transcript expression levels 
of OsCaml-1, and especially of OsCML4, 5, 8, and 11 
were significantly increased under dehydration, whilst 
that of OsCML 13 was slightly down-regulated. Under 
salt stress, significantly increased transcript expression 
levels of OsCarnl-1, and especially of OsCML4, 8, 9, and 
11 were observed, whilst that for OsCML13 was slightly 
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Table 1 Oligonucleotide primers used in this study 



Primer name 


Sequence 


Position 3 


Amplicon size (bp) b 


Annealing temp. (°C) 


OsCam 7- 7-F 


5'- ACCGTGCA^GCCGTATOG -3' 


499-518 


177 


58.3 


OsCam 1- 7-R 


5'- GCAAGCCTOACAGATOAC -3' 


656-675 






OsCML 7-F 


5'- CCAGAAGTGCGTGATCCTGT -3' 


543-562 


184 


58.3 


OsCM 7-R 


5'- ACTACGGACTACGGCTGTGA -3' 


707-726 






OsCML3-F 


5'- ACTACAACGAGTOCTCAAG -3' 


410-429 


180 


57.3 


OsCML3-R 


5'- CATCAGAACAGTOCAAACC -3' 


570-589 






OsCML4-F 


5'- GCAGGTGAACTACGATGAAT -3' 


402-421 


193 


56.3 


OsCML4-R 


5'- TACCCATAGCTGAAGTCCAA -3' 


575-594 






OsCML5-F 


5'- ATGATGCTCTCCGACCAATA -3 ; 


481-500 


180 


57.3 


OsCML5-R 


5'- CCAAGGCCAAATOAATGAC -3 ; 


641-660 






OsCML7-F 


5'- CCGCATCGTCGCCAAATAAT -3' 


429-448 


193 


57.3 


OsCML7-R 


5'- GTCCAAATCACACCGGAATG -3' 


602-621 






OsCML8-F 


5'- AGATGATGAAGAGGATAGGA -3' 


539-555 


185 


56.3 


OsCML8-R 


5'- AAACATAAGGCGGTATGGTA -3' 


701-720 






OsCML9-F 


5'- TACAAGGAGTOGTCAAGGT -3' 


430-449 


170 


58.3 


OsCML9-R 


5'- GATOGOTGAATCATATCGC -3 ; 


579-600 






OsCMLl 1-V 


5'- CAACATOTCTCCTGAGAAT -3' 


621-640 


183 


56.3 


OsCMLll-H 


5'- ATOACAAGAGCTCGATCAC -3' 


784-803 






OsCMLl 3-V 


5'- ATCGAAATGGTGATGGTGAG -3 ; 


437-456 


193 


58.3 


OsCMLl 3-R 


5'- GCATGGTOTOTOTOAG -3' 


610-629 






OsEFla-F 


5'- ATGGTOTGGAGACCTO -3' 


1192-1209 


127 


58.3 


OsEFla-R 


5'- TCACOTGGCACCGGTO -3' 


1301-1318 







a Position of the primers from the GenBank sequence given in the materials and methods, where position 1 is the predicted open reading frame start codon and 
numbered 5' to 3' on the sense strand. 

b Expected amplicon size based upon the primer positions on the GenBank sequence. 



down-regulated. Note that although numerically the 
transcript expression level of OsCML3 was increased fol- 
lowing both dehydration and salt stress, this was not 
statistically significant due to the large variance in the 
control samples, and so remains to be clarified. 

With respect to the three-week-old O.sativa 'KDML105' 
seedlings, rt-RT-PCR revealed that under osmotic stress 
(20% PEG), OsCML4 and OsCML8 exhibited a nine- and 
four- fold increases in their transcript expression levels, 
respectively, and that this started from 3 h after osmotic 
stress and were maintained at these elevated expression 
levels until at least 48 h after treatment (Figure 3B). How- 
ever, whilst significantly elevated transcript expression 
levels of OsCML 11 were also detected at 3 h after os- 
motic stress (13-fold), in contrast these expression le- 
vels then declined to six-fold higher and control levels 
at 6 h and 12 h onwards after osmotic stress, respecti- 
vely. Expression of OsCaml-1, and OsCMLS showed an 
earlier response, being up-regulated 5-6-fold at 1 h after 
osmotic stress treatment compared with the control, but 
whilst OsCaml-1 expression levels remained slightly up- 
regulated at 3 to 48 h after osmotic stress that for 
OsCMLS returned to their normal levels after 3 to 6 h. 



OsCML7 and OsCML9 exhibited a slight increase in their 
transcript expression level at early time points (at 1-3 and 
3-6 h after treatment, respectively), while no significant 
changes in the transcript expression levels of OsCMLl, 3 
and 13 were observed at all time points examined. 

Similar to osmotic stress, OsCML4 and OsCML8 also 
showed significant prolonged increases in their trans- 
cript expression levels under salt stress (150 mM NaCl), 
starting from 3 h until at least 24 h after treatment 
(Figure 3C). Compared to the respective control, the 
OsCML4 transcript expression levels were maximal at a 
four-fold elevated level at 6 and 12 h after treatment 
while OsCML8 reached the highest levels (12-fold) at 12 
and 24 h after treatment. Early up-regulated transcript 
expression of OsCaml-1 and OsCMLS at 1 h after treat- 
ment (five- and eight-fold, respectively) and of OsCMLl 1 
at 3 h after treatment (six-fold) was also observed. In 
contrast, the transcript expression levels of OsCML3 
and OsCML13 were not significantly altered in res- 
ponse to salt stress at all time points examined in this 
study, whilst OsCML7 and OsCML9 were only slightly 
and transiently increased at 1 h and 6 to 12 h after salt 
stress, respectively. However, unlike under osmotic stress, 
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Figure 2 OsCc/m and OsCML transcript expression levels in various 'KDML105' rice organs as determined by rt-RT-PCR. Relative transcript 
expression levels of OsCaml-1 and the nine OsCML genes in (A) leaf blades/sheaths of three-week old 'KDML105' rice seedlings relative to 
0sCML5 and in (B) leaf blades/sheaths, roots, flowers, and seeds of the 'KDML105' rice during the grain filling period compared with their 
respective level in the leaf of the three-week old rice seedlings. Data are shown as the mean ± 1 SD, and are derived from three independent 
replicates. Means with a different lowercase letter are significantly different (p < 0.05). 



OsCMLl exhibited increase in its expression at 1 h after 
salt stress. 

Examination of the OsCaml-1 and OsCML promoters 

In order to examine the promoters of these genes and 
their ds-acting elements, the sequences upstream of 
their coding regions were retrieved from GenBank and 
bioinformatically analyzed using the PLACE [13,14] and 
the PlantPAN [15] software and dtabases. Several cha- 
racteristic elements were located; including motifs invol- 
ved in responses to dehydration, abscisic acid (ABA), salt 
stress, and cold stress (Figure 4). Putative DREs (drought- 
responsive elements) [16] were found in the 5' flanking 
(upstream) region of five of the loci (OsCaml-1 at -1064, 



OsCML4 at -393, OsCMLS at-892, OsCML7 at-954, and 
OsCMLl 1 at-561 and -565). In addition, the 5' upstream 
sequences from eight of these genes (OsCaml-1, OsCMLl, 
3, 4, 5, 7, 8 and 13) contained putative ABREs (ABA- 
responsive elements) [17], while LTREs (low temperature 
responsive elements) [18] were found in seven of these 
genes (OsCaml-1, OsCMLl, 3, 4, 5, 7, and 11). 

To investigate inducibility of the OsCaml-1 promoter 
by salt-stress, (0 (control), 100, 150 or 300 mM NaCl), the 
(3-glucuronidase (GUS) activity of three independent lines 
of the homozygous T3 transgenic KDML105 rice cultivars 
harbouring the OsCaml-lr.gus construct was examined 
fluorometrically in two-week-old seedlings. All three inde- 
pendent transgenic lines exhibited similar trends of GUS 
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Figure 3 (See legend on next page.) 
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f N 

(See figure on previous page.) 

Figure 3 OsCam and OsCML transcript expression in response to stress. (A) Microarray based transcript expression profiles under 
dehydration (drought) and salt stress shown as a heat map. Transcript expression levels examined in the KDML105 rice cultivar under osmotic 
stress (B) and salt stress (C) showing relative transcript expression levels standardized to that of OsEFla and expressed relative to the levels on 
day 0 of the treatment. Rice seedlings were grown in a CRD with three replicates, and for each replicate, ten seedlings were pooled for RNA 
extraction. The PCR reaction of the same cDNA preparation was performed in triplicate for technical replication. Data are shown as the mean ± 1 
SD and means with a different lowercase letter are significantly different (p < 0.05). 



activity to each other under the different NaCl treatments 
(Figure 5). The plants treated with 100 mM NaCl did not 
exhibit any statistically significant difference in GUS acti- 
vity in all three independent lines (p < 0.05) compared to 
their respective controls at 0 and the untreated plants, 
whereas at 150 mM NaCl all transgenic lines had exhi- 
bited a significantly higher GUS activity from 4 h after 
treatment in lines 1 and 3, and after 1 to 8 h in line 2. For 
the plants treated with 300 mM NaCl, a biphasic induc- 
tion of GUS activity was observed in all three independent 
transgenic lines with an initial induction and highest in- 
duction of GUS activity at 1 h after treatment, which then 
decreased to the control level by 2 h and was followed by 
the second phase of induction at 4 h after treatment. This 
second phase was either transient (line 2) or maintained 
for up to 8 h (line 1) or 24 h (line 3) after treatment. 

Discussion 

Whilst the OsCaM group contains five members that 
share a highest degree of amino acid sequence identity 
(> 97%) to known typical CaMs from other plants, pre- 
vious report phylogenetically classified the highly con- 
served OsCML proteins into four groups [5]. The OsCML 
proteins are small proteins of 145 to 250 amino acid 



residues with an approximately 44% to 85% amino 
acid identity to typical plant CaMs. A large family of 
six Cam and 50 CML genes have also been annotated 
in the A.thaliana genome [8], suggesting the existence 
of an extensive set of CaM and CML proteins in each 
plant species. 

Even though rice is considered an important crop and 
a model for other monocots and especially cereal crops, 
almost all of OsCam and OsCML genes have not been 
characterized. From the publicly available microarray 
data, all OsCam genes were found highly expressed in 
almost all organs examined (Figure 1A). One of the de- 
fining characteristics of CaMs in plants is the presence 
of multiple CaM isoforms. Even though they are ubiqui- 
tously expressed, the different isoforms can display dif- 
ferential responses to individual stimuli in any given 
plant tissue [5,11,19-22], suggesting that each of the 
OsCam genes may have distinct physiological functions 
depending on where and how the expression of each 
gene is up-regulated in addition to the different bio- 
chemical properties that might be affected by the slight 
differences in their primary structures [21,23]. 

In contrast, several highly conserved OsCML genes 
(groups 2-5) genes tend to be highly expressed in some 
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organs/tissues and their levels of expression are modu- 
lated during different stages of development (Figure 1A), 
indicating that they are developmentally regulated. It is 
conceivable that CaM proteins have many targets and 
are ubiquitously involved in numerous cellular processes 
while the highly conserved CML proteins have more 
specialized targets and are involved in more specific pro- 
cesses of the cell. Interestingly, each subgroup of the 
highly diverged group of OsCML (group 6a to 6e) genes 
had one member whose expression is found at a notice- 
ably higher level than the other members. Members in 
the same subgroups have the same number and config- 
uration of EF-hand motifs. It is speculated that these 
highly expressed genes may be the predominant gene 
representing each subgroup that has essential functions 
in the cell. 

When the microarray data was closely examined for 
the transcript expression levels of the OsCaml-1 and 
nine OsCML genes (Figure IB), ubiquitous expression at 
relatively high levels of OsCarnl-1, OsCML7, and 13 was 
observed, which suggests that they may have important 
functions during the regular growth and development 
processes of rice plants. The increasing transcript ex- 
pression levels of OsCMLl, 4, and 8 with maturation of 
the inflorescence until they reached similar levels as 
those in the leaf and root tissues in the early developing 
seed, suggest that they may have important roles during 



the maturation of the inflorescence and the early seed 
development, especially OsCMLl, which expression level 
peaked during those stages. 

The transcript expression levels in the 'KDMLIOS' rice 
in the leaf blades/sheaths of the three-week old 
'KDML105' rice seedlings, as evaluated by rt-RT-PCR 
(Figure 2A) agree well with the microarray data in the 
IR64 cultivar (Figure 1A). When the transcript expres- 
sion levels of OsCaml-1 and OsCMLs were examined in 
different organs (Figure 2B), their profiles could be 
divided into three groups as those with the highest ex- 
pression level in (i) the root, or (ii) in the flower and the 
seed, or (iii) those genes with similar expression levels 
among the organs examined. Relatively high transcript 
levels of OsCML3 and OsCMLll in the flower and the 
seed, and OsCML4 } and 5 in the root suggest their 
functional significance in those respective organs. The 
OsCaml-1 and the other five OsCML genes {OsCMLl, 
7, 8, 9 and 13) had more or less similar expression levels 
among the different organs examined. From the micro- 
array data of the IR64 rice cultivar, four of these genes 
(OsCarnl-1, OsCML7, 9 and 13), were expressed at more 
or less constant levels in different developmental stages. 
However, expression patterns among different organs/ 
tissues that are not consistent with those analyzed from 
the microarray database may also result from variation 
in the rice varieties examined. 
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Up-regulated expression of a gene in response to a 
stress signal may reflect the function of the correspon- 
ding gene product, especially in signal cascades. In a 
large gene family, investigating expression patterns of 
their members could point to genes or isoforms that po- 
tentially function under the conditions of interest. In this 
study, the transcript expression analysis of OsCaml-1 
and nine OsCML genes in the KDML105 rice cultivar by 
rt-RT-PCR revealed that expression levels of OsCaml-1, 
OsCML4, 5, 8, and 11 were increased under osmotic 
stress (20% (w/v) PEG) (Figure 3B), which is consistent 
with the GSE6901 data for the IR64 cultivar from the 
DNA microarray database (Figure 3 A) [12] and suggests 
that these genes may function in the mechanisms of 
Ca 2+ -mediated responses to osmotic stress. Similarly, ex- 
pression of these genes was also found to increase in 
the KDML105 rice cultivar under salt stress (150 mM 
NaCl) along with OsCMLl (Figure 3C), suggesting that 
OsCaml-1, OsCMLl, 4, 5, 8, and 11 may function in the 
mechanisms of Ca 2+ -mediated responses to salt stress. 
The up-regulated transcript expression levels under salt 
stress of OsCML4, 8 and 11 were also consistent with 
the GSE6901 data for the IR64 cultivar from the DNA 
microarray database (Figure 3A). Whether the discrep- 
ancy in the transcript expression patterns of OsCaml-1, 
OsCMLl, and 5 is due to the different timings of expres- 
sion being monitored or differences in the rice varieties 
examined remains to be evaluated but, overall, up- 
regulation of almost the same set of the highly con- 
served OsCML genes under osmotic and salt stresses 
was observed, confirming that conclusion that genes re- 
sponsive to osmotic stress overlap to a high extent with 
those that are responsive to salt stress [24]. 

Several of these genes not only exhibited up-regulated 
transcript expression levels by osmotic stress and salt 
stress, but also exhibited different patterns of up- 
regulation in terms of timing and levels of expression 
within a tissue type and exhibited differential expression 
in different tissues/organs. Together, differential temporal 
and spatial expression patterns of these OsCaml-1 and 
OsCML genes suggest that each individual gene product 
may possess specific roles during Ca 2+ -mediated responses 
to osmotic stress and salt stress. Comparison of the tran- 
script expression patterns of these genes under osmotic or 
salt stress in the KDML105 rice cultivar, as determined by 
rt-RT-PCR in this study, with the publicly available RNA- 
Seq data and the other microarray data sets [25], revealed 
that OsCML4 and OsCML8 consistently exhibited higher 
expression levels under osmotic and salt stresses suggest- 
ing their significant functions in Ca 2+ -mediated responses 
to these stimuli. In addition, the early induction within 1 h 
of OsCaml-1 and OsCMLS under osmotic and salt 
stresses, and of OsCMLl under salt stress suggests 
their importance in conveying the stress signals early in 



the transduction cascades of Ca + signaling. However, cau- 
tion must be taken in interpreting these results because 
the changes in levels of transcripts generally do not co- 
incide with the changes in the levels of the proteins they 
encode and, given the time it takes to express a protein, 
may not necessarily reflect their involvement in the 
responses of the cell to stress. Nonetheless, a gene of 
which transcript expression is early induced would be 
given special attention as a candidate for further inves- 
tigation into its possible involvement in response to a 
particular stress. 

Several cis-acting elements in response to abiotic 
stresses in the 5' upstream regions of the OsCaml-1 and 
OsCML genes were detected (Figure 4). DREs, important 
promoter elements that are responsive to drought, high 
salt, and cold [16], were located in the putative promo- 
ters of OsCarnl-1, OsCML4, 5, 7, and 11. Consistent 
with this is that the transcript expression level of all 
these genes was shown by rt-RT-PCR to increase in the 
KDML105 rice cultivar under osmotic and salt stresses 
with OsCML7 exhibiting a slight increase in its expres- 
sion level and the others significantly much higher level 
of up-regulation. DREs specifically interact with the 
transcription factors DREBs and regulate expression of 
many stress-inducible genes. In rice, the DRE binding 
protein 1 (OsDREBl) functions in the cold stress res- 
ponse, whereas OsDREB2 functions in the heat and os- 
motic stress responses [26]. The results here indicate 
that these DREs are potentially responsible for the os- 
motic stress-induced expression of these genes and 
interesting candidates for further characterization. 

In the 5' upstream region of the OsCaml-1 promo- 
ter, the putative DRE motif (ACCGAC) was located at 
-1062. Induction of the GUS activity level in the three 
independent OsCamwgus transgenic rice plants (Figure 5) 
suggested that induction of OsCaml-1 expression under 
salt stress is, at least, partly due to the activity of its pro- 
moter. Its induction appears to be dependent on the se- 
verity of the salt stress (concentration of NaCl), with no 
significant induction in the plants treated with 100 mM 
NaCl while a monophasic and biphasic induction were 
observed in the plants treated with 150 mM and 300 
mM NaCl, respectively, in all three independent trans- 
genic lines. A biphasic induction of OsCaml-1 expres- 
sion by salt stress in the KDML105 rice cultivar, as 
determined by rt-RT-PCR has been reported previously 
[27], whilst heat shock (HS) induced biphasic [Ca 2+ ] cyt 
signal in rice root cells and the HS -induced expression 
of OsCaml-1 strongly oscillated [11]. The complex 
responses of the OsCaml-1 gene to salt stress suggest 
that OsCaml-1 is a significant player in the Ca 2+ signal 
transduction network under salt stress. However, cau- 
tion must be taken in interpreting these results, espe- 
cially on timings of the induction because of the possible 
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differences in mRNA and protein stability between the gus 
gene and the OsCaml-1 gene. 

Conclusions 

Differential expression patterns in tissues/organs were re- 
ported among different Cam and CML genes in rice. Un- 
der osmotic and salt stresses, expression of several genes 
including OsCaml-1, OsCML4, 5, 8, and 11 was induced 
with different kinetics and magnitude. In agreement with 
their inducibility, several ds-acting elements in response 
to abiotic stresses including DREs, important promoter 
elements responsive to drought, high salt, and cold, were 
detected in the 5' upstream regions of these genes. Induc- 
tion in the GUS activity of the transgenic rice expressing 
gus gene under the control of the OsCaml-1 promoter, 
which contains a putative DRE motif (ACCGAC) at 
-1062 was observed. Its induction appeared to be biphasic 
and dependent on the severity of salt stress. These results 
suggest that these OsCam and OsCML genes play diffe- 
rential roles as sensor relays in regulating developmental 
processes and Ca 2+ -mediated responses to abiotic stress. 

Methods 

Materials 

TRI REAGENT was purchased from Molecular Research 
Center, Inc. (Cincinnati, OH, USA). The iScript cDNA 
synthesis kit and the SsoFast EvaGreen Supermix were 
purchased from Bio-Rad (Hercules, CA, USA). 4- 
methylumbelliferyl p-D-glucuronide (4-MUG) and 4- 
methylumbelliferone (4-MU) were from Sigma (St. 
Louis, MO, USA). Synthetic oligonucleotides for real- 
time reverse transcription polymerase chain reaction 
(rt-RT-PCR) were obtained from 1 st Base (Singapore). 
Seeds of Oryza sativa L. cultivar Khao Dawk Mali 105 
(KDML105) were provided by the Rice Research Center 
(Patumthani, Thailand). 

Plant growing and stress treatments 

O. sativa 'KDML105' seedlings were hydroponically grown 
in nutrient solution under a 12-h light/ 12-h dark photo- 
period. After three weeks, the plants were treated with 
150 mM NaCl for 1, 3, 6, 12, 24 or 48 h by adding the 
same medium except supplemented with 150 mM NaCl 
or 20% (w/v) PEG4000 into the chamber containing the 
rice seedlings with minimal touch contact and distur- 
bance. Rice seedlings were grown with a completely ran- 
domized design (CRD) with three replicates, and for each 
replicate, ten seedlings were pooled for RNA extraction. 

Real-time reverse transcription polymerase chain reaction 
(rt-RT-PCR) 

Leaf blades/sheaths of three-week old seedlings; or leaf 
blades/sheaths, roots, flowers and seeds during the grain 
filling period were collected and immediately frozen in 



liquid nitrogen and stored at -80°C. Rice tissues were 
ground in liquid nitrogen to a fine powder using a 
chilled mortar and pestle. Total RNA from the frozen 
tissues was extracted using the TRI REAGENT with 0.2 
mL chloroform and precipitated by mixing with isopro- 
panol. Subsequently, cDNA was synthesized by RT-PCR 
using the iScript cDNA synthesis kit according to the 
manufacturer s instruction. Then, rt-PCR was performed 
in a final volume of 20 |iL, which contained a 2-(iL ali- 
quot of the first strand cDNA reaction, 0.05 mM of each 
of the gene-specific primers (Table 1), and lx SsoFast 
EvaGreen Supermix. The reaction included an initial 8 
min at 95°C, followed by 40 cycles of 95°C for 30 s; 
Ta°C for 30 s and 72°C for 45 s (where the annealing 
temperature, Ta, was as in Table 1). The specific oligo- 
nucleotide primers for the second stage of the rt-RT- 
PCR (rt-PCR of the cDNA), shown in Table 1, were 
selected using the Primer3 algorithm [28] from the 
OsCaml-1 and OsCML cDNA sequences (see sequence 
analysis below). 

The O.sativa elongation factor- la (OsEFla) was as- 
sumed to be independent of the abiotic stresses in its ex- 
pression levels (a housekeeping gene), and so amplified as 
an internal control using the specific primers OsEFl<x-F/R 
(Table 1), designed as above against the corresponding se- 
quence (GenBank accession code AK105030). The level 
of transcript expression of the OsCaml-1 and respective 
OsCML genes expression was then standardized against 
that for the OsEFla gene transcript expression level from 
the same cDNA template as an internal control, and then 
expressed with reference to the standardized transcript ex- 
pression level on day 0 of the treatment. Rt-RT-PCR reac- 
tions of each cDNA preparation were performed in 
tripicate for technical replication. Data were compared 
using the analysis of variance (ANOVA), and then the 
means were compared with Duncans multiple range test 
(DMRT) accepting significance at the p < 0.05 level. 

Determination of (^-glucuronidase (GUS) activity 

The homozygous seeds of three independent OsCaml-1: 
gus transgenic rice lines [10] were hydroponically grown 
in 0.25 x strength NB medium [29] under a 16-h light/ 
8-h dark photoperiod and treated with salt stress by add- 
ing, with minimal contact and disturbance, the medium 
containing 0 (control), 100, 150 or 300 mM NaCl into 
the chamber containing the rice seedlings to be treated. 
The leaf tissues of two-week-old rice seedlings were col- 
lected and frozen in liquid nitrogen at 0, 1, 2, 4, 8, 24 
and 48 h after treatment. Rice seedlings were grown in a 
CRD with four replicates, and for each replicate, two 
seedlings were pooled for RNA extraction. The collected 
tissues were ground with a chilled mortar and a pestle, 
and homogenized in protein extraction buffer (50 mM 
Sodium phosphate, pH 7.0, 10 mM EDTA, 10 mM |3- 
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Mercaptoethanol, 0.1% Sodium n-lauroylsarcosine, 0.1% 
Triton X-100). The p-glucuronidase enzyme (GUS) reac- 
tion assay was performed in protein extraction buffer 
containing 1 mM 4-MUG. For each reaction the sub- 
strate solution (200 \iL) was pre-incubated at 37°C for 5 
min, and then the reaction was initiated by the addition 
of 1 \iL of the test plant extract and incubated at 37°C. 
After 10, 20, 30 and 40 min of incubation time 10 \iL ali- 
quots of the reaction mixture was transferred to 100 \iL 
of stop solution (0.2 M Na 2 C0 3 ). The standard curve 
was generated using five concentrations (10 nM to 100 
nM) of 4-MU. The fluorescence intensity of the samples 
was measured by setting the excitation at 365 nm and 
the emission at 455 nm. The p-glucuronidase activity of 
the protein extracts was calculated as pmole 4-MU pro- 
duced per minute per \ig protein (units per \ig protein). 
Fluorescence of each replicate was measured three times. 
Protein concentrations of the samples were determined by 
the Bradford's method [30] using bovine serum albumin 
as a standard. Data were compared using ANOVA, and 
the means were compared with REGWQ accepting signifi- 
cance at the p < 0.05 level. 

Sequence analysis and microarray data retrieval 

To determine which OsCML genes may have altered ex- 
pression patterns in response to stress, their expression 
levels in different rice organs and tissues were compared. 
DNA microarray data of the OsCam and OsCML genes 
were retrieved from the publicly available microarray 
data of global gene expression from the experiment GSE 
6893, which examined rice gene expression in various 
rice tissues/organs and stages of reproductive develop- 
ment from the IR64 rice cultivar or GSE 6901, which 
examined gene expression of 7-day-old 1R64' rice seed- 
lings grown under 3 h of dehydration (drought) or salt 
stress (200 mM NaCl) [12] via the rice oligonucleotide ar- 
ray database [25] available at (http://www.ricearray.org/). 

For PCR design, nucleotide sequences from O.sativa 
were retrieved from of the Rice Genome Annotation 
Project database [31] and GenBank at the National 
Center for Biotechnology Information (http://www. 
ncbi.nlm.nih/). Sequences of the following genes (loci): 
OsCaml-1 (LOC_Os03g20370), OsCMLl (LOC_Os01g5 
9530), OsCML3 (LOC_Osl2g03816), OsCML4 (LOC_Os 
03g53200), OsCMLS (LOC_Osl2g41110), OsCML7 (LOC_ 
Os08g02420), OsCML8 (LOC_Osl0g25010), OsCML9 
(LOC_Os05g41200), OsCMLll (LOC_Os01g32120) and 
OsCML13 (LOC_Os07g42660) were obtained. 

To identify ds-acting regulatory DNA elements with- 
in the promoters of the OsCaml-1 and OsCML ge- 
nes, their 5' upstream sequences were analyzed using a 
Database of Plant Cis-acting Regulatory DNA Elements 
(PLACE) (http://www.dna.affrc.go.jp/PLACE/) and the 



Plant Promoter Analysis Navigator (PlantPAN) (http:// 
plantpan.mbc.nctu.edu.tw/) . 
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